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Abstract 

We use the detection of neutrinos from a future galactic type-II supernova event in a water 
Cerenkov detector like Super-Kamiokande to constrain the possible violation of spin-statistics 
by neutrinos resulting in their obeying a mixed statistics instead of Fermi-Dirac. 
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Recently there was a suggestion of neutrinos violating the spin-statistics relation and thereby 
becoming a good candidate for all (or part) of the dark matter in the universe [1]. There are a 
number of papers which discuss the possible violation of spin-statistics relation by neutrino though 
no consistent satisfactory model exists [2, 3, 4, 5]. On the other hand the experimental verification 
of neutrinos with half spin following a mixed statistics was not thoroughly studied earlier. For 
nucleons however such studies exist [2]. The recent paper on the possibility of neutrinos violating 
the Pauli Exclusion Principle [1] however has renewed interest in the subject [6]. The Doglov- 
Smirnov work [1] points out that if neutrinos do obey Bose-Einstein (BE) statistics instead of 
Fcrmi-Dirac (FD), then they may form large cosmological Bose condensates and account for the 
dark matter. This also opens up the possibility of large lepton asymmetry in the universe. As 
double beta decay disallows purely bosonic neutrinos, one will be interested in mixed statistics in 
terms of a continuous "Fermi-Bose" parameter, k {k — —1 is purely fermionic and k = 1 is purely 
bosonic) [7]. This on the other hand has astrophysical consequences as well. For example this 
will have impact on the type II supernova (SN) dynamics and will change the energy spectrum 
of neutrinos coming out of the supernova. In this report we shall not be concerned with the 
justification of the Dolgov-Smirnov suggestion but concentrate on the question: Given the scenario 
of neutrinos of all three flavors obeying the mixed statistics, how would future observation of 
galactic neutrinos in large terrestrial detectors put limits on the mixed statistics parameter k. We 
assume that apart from (possibly) violating the spin-statistics theorem, the massive neutrinos do 
not have any other non-standard property. We shall see that with detectors like Super-Kamiokande 
(SK) one indeed can test this hypothesis at a significant level of accuracy. 

Massive stars at the end of their normal lifespan, collapse due to the gravitational pull once 
the silicon burning in the core stops and the core has a mass greater than the Chandrasekhar 
mass. The neutrinos that are produced by electron capture on nuclei and free protons during the 
initial collapse phase escape. However, as the density of the core exceeds densities of 10^^ — 10^^ 
gm/cc during collapse, neutrinos get trapped. At densities higher than the nuclear matter density 
a shock wave forms inside the core and travels outward. Whether the shock wave can reach the 
edge of the core with enough energy to cause the explosion with the observed energies is the central 
question of supernova physics today. 

At the high densities and temperatures of the SN core during the post-bounce phase, neutrinos 
and antineutrinos of all three flavors are produced. Almost ~ 99% of the gravitational energy 
produced through the huge contraction of the star (~ a few times 10^^ ergs) is released in the form 
of the six species of neutrinos {ve, J^n, and their antiparticles) . In the high density central part 
of the core, ncTitrinos of each flavor have high opacity and hence cannot come out. But at larger 
radius due to a strong density gradient, neutrinos can diflusc out and eventually escape. Therefore, 
we expect almost thermal spectra for all the neutrino species, with temperatures characteristic of 
their radius of last scattering in the SN, usually called the neutrino sphere radius. Since and Pg 
have charged current interactions with the SN matter in addition to neutral current interactions 
and since SN matter is neutron rich, z/g, Ve and z/^ ('x' stands for r, // and f ) decouple at different 
radii. Hence the neutrino spheres of the three different types, z/gj and Vxi have different radii 
and hence different equilibrium temperatures (Tj,), with T^,^ < Tp^ < T^^{Tx). For our purpose 
we shall take T^^ = 3.5 MeV, Tp^ = 5 MeV and T^^ = 8 MeV [8], used in most simulations 
and calculations as the standard values. As there are arguments which claim that Tp^ and T^^ 
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should be closer [9], we vary these temperatures over a range as detailed later in this report. 
These temperatures are obtained by using FD equilibrium energy distributions and may change 
somewhat had the simulations been carried out with BE distributions or distributions with mixed 
statistics. However the qualitative results of our investigation here will not change if they are done 
with equilibrium BE temperatures for z/g, and v^. Again realistic simulations indicate small 
departures from equilibrium distribution of energies of the neutrino, with the high energy tail 
lower or "pinched" . The pinching factor^ is expected to range between 77e~0 — 3, r/g^O — 3 and 

~ — 2 for the fg, and flux respectively [10]. We further assume that the total luminosity 
Eb is equipartitioned among the six neutrino and antineutrino species. We remind readers here 
that the electron type antineutrinos (with at most one neutrino scattering event) were detected 
alongwith the SN 1987A explosion. But with 11 and 8 events at Kamioka and 1MB, fitting of 
energy spectrum and determination of temperature had large errors [11]. 

Following [1] we parametrize the equilibrium distribution of the neutrinos with a mixed statis- 
tics as 



where E is the energy of the neutrinos, T the temperature of the distribution and k is the Fermi- 
Bose parameter for the mixed statistics, as mentioned earlier. Fig. 1 shows the energy spectrum of 
the neutrinos corresponding to FD [k = +1), Maxwell Boltzmann MB [k = 0) and BE {k = — 1) 
distributions, with typical temperatures of 3.5 (left panel), 5 (middle panel) and 8 MeV (right 
panel) . The area under each of the curves gives the average energy expected for the corresponding 
values of T,^ and k. It is clear from the figure that for the same T^, at small values of E the 
predicted spectrum for the BE distribution is higher than for the FD distribution. However for 
large E, the trend is reversed and the spectrum predicted for FD is higher than that for the BE 
distribution and has a longer high energy tail. This trend reflects the fact that for the same 71^, 
the predicted average energy for the neutrinos is larger for the FD distribution. We also stress 
from Fig. 1 the fact that the spectral shape for the FD and BE distributions arc different. We 
want to use the signature of this difference in predicted average energy and spectral shape of the 
SN neutrinos in terrestrial detectors, to disfavor the "wrong" neutrino statistical distribution and 
put limits on the Fermi-Bose parameter k. 

In Fig. 2 we show the effect of introducing the pinching in the FD distribution. We plot the 
galactic SN neutrino flux spectrum expected at earth given by [10] 



where is the SN luminosity released in D is the distance of the SN (10 kpc), r]a is the 
pinching factor and 



fjf"^ = (exp(E/T) + k) 



(1) 




La E'^ 



(2) 




(3) 



^The definition of t] is given in Eq. (2). 



3 



The black solid line shows the spectrum for a pure FD distribution with rja = 0, the red dotted line 
corresponds to a pinched FD distribution with rja — "2 while the green dot-dashed hne is for the 
pure BE distribution. All the three spectra plotted in the figure have the same average energy^. 
A comparison of the three spectra shows that the difference in the spectral shape between the 
BE and FD distribution increases as rja changes from zero. Hence, we expect it would be easiest 
to confuse between the BE and FD distributions when r^a = and it is for this case that we 
generally expect the minimum x^- Since we want to check how well one can differentiate between 
the FD and BE (and mixed) distributions and since this seems to be most difficult for the pure 
FD distribution, we will neglect the effect of pinching in the spectrum in the rest of the paper for 
simplicity'^. 

We shall be concerned with the large water Cerenkov detectors like SK. The main interaction 
channel in SK which is important for the detection of the SN neutrinos is the charged current 
capture of Pe on free protons: 

Pe + P — e"*" + n (4) 

with a threshold 9e energy of 1.3 MeV. The number of positrons expected in SK from a galactic 
SN explosion is given by 

R^NjdEAj dEr R{Ea, Et) J dE ap^p{E) [Fe{E) P,-,{E) + F^{E){1 - Pee-(^))] , (5) 

where N is the number of target protons and (Tu^p is the z/g-p capture cross-sections. E is the 
energy of the incident Ve-, while Ea and Et are respectively the measured and true energy of the 
emitted positron. The true and measured energy of the positron are related through R{Ea, Et), 
the Gaussian energy resolution function of the detector given by. 



R{Ea,Et) = -^=exp \ , (6) 




where we assume ctq to have a value similar to what is used in the analysis of the SK solar neutrino 
data. Eg and F^ are the flux of and produced in the SN core and depend on the spin statistics 
of the neutrinos and Pge is the survival probability of the SN z/g. 

Before reaching the detector, the neutrinos travel through the SN matter, propagate through 
vacuum and finally may even travel through earth matter depending on the time of the SN burst. 
For the range of neutrino oscillation parameters consistent with the solar and KamLAND reactor 
data (Am| = Aml^ = 8.0 x 10"^ eV^ and sin^ 9q = sin^ 9^ = 0.31 [12, 13]) and SK atmospheric 
data (Am^tj^ = Aml^ = 2.1 x 10"^ eV^ and sin^ 2^^atm = sin^ 2^^23 = 1-0 [14]), one expects large 
matter enhanced oscillation inside the supernova and possible regeneration effects inside the earth 
[10, 15]. As a result, the more energetic and P,- get converted to Pe and make the spectra harder. 
This results in enhancing the number of expected charged current events in the detector and is 

^Note that this is different from Fig. 1 where the temperature was kept same for the different distributions. 

^Thc issue of pinching could get compUcated in some cases, especially for a BE spectrum with chemical potential 
and in these cases one has to be very careful. We propose to make a detailed study of this in a future work. However, 
we have checked that taking no pinching in the spectra yields the smallest for most realistic cases. 
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incorporated in our numerical analysis through Pgg in Eq. (5). We assume that the supernova 
explodes at a distance of 10 kpc from us and that the time of the event is such that the neutrinos 
do not cross the earth matter. The survival probability Pge depends critically on the value of ^13, 
which at the moment has a very weak bound of sin^ ^13 < 0.04 at 3a [12, 16]. In what follows, we 
will assume that the true values of the oscillation parameters Arri^i, sin^ 612 and Am|]^ have been 
determined with sufficient accuracy to the values mentioned above (see [17] and [18] for recent 
detailed discussions). For the major part of the paper, we will also assume that the true value of 
mixing angle ^13 ~ 0. We will discuss the impact of the uncertainty on this oscillation parameter 
towards the end. 

Fig. 3 displays the distribution of the number of (z/g + p) events in SK as a function of the 
measured energy Ea of the positron. The black solid (dashed) line gives the event spectrum for 
a FD distribution with (without) neutrinos oscillations, while the red sohd (dashed) line gives 
the expected event spectrum for a BE distribution. We assume Tp^ = 5 MeV, — S MeV and 
Eb = 3 X 10^^ ergs. From the figure, BE distribution can be seen to undcrprcdict the number of 
events over a large range of energy from about 12-70 MeV, compared to FD distribution. We also 
give the ilcr error bars corresponding to the number of events SK would observe in each bin for 
the FD distribution. This clearly shows that a large detector like SK should be able to distinguish 
between the different statistics at a good confidence level. 

In order to quantify our statement, we perform a statistical analysis of the neutrino induced 
positron spectrum that SK would observe in the event of a galactic SN explosion. For the errors 
we take a Gaussian distribution and define our as 

= Y^iNf'^'^" - Nf'^'^){al)-\Nf"^y - Nf'^'^) , (7) 

where _/V*'*^°^^ ^jyobsvd^ ^^le number of predicted (observed) events in the i^^ energy bin and 
the sum is over all bins. The error matrix afj contains the statistical and systematic errors. The 
systematic uncertainty will be eventually determined by the experimental collaboration after a 
real SN event has been observed in SK. However, for the moment we have simply assumed a 5%^ 
systematic uncertainty in the simulated "data", fully correlated in all the bins. The number of 
energy bins used in our analysis is shown in Fig. 3. The function is minimized with respect 
to the parameters in the theory to give us a handle on the CL. with which observation of SN 
neutrinos in SK could lead to the determination of their correct distribution function. In what 
follows, there are two kinds of statistical tests that we shall perform. In one, we shall compare the 
data set obtained for (say) the true FD distribution with the predicted event spectrum for the false 
BE distribution. This will give us the measure of the C.L. with which the false BE distribution 
could be ruled out. In another, we will compare the data set generated for (say) the true FD 
distribution (/t = 1) and compare it with a distribution for all possible values of k in the range 
[— 1, 1], with K = —1 corresponding to a pure BE distribution. 

In Fig. 4 we show the x^ ^ ^ function of the parameter k. We assume that the only errors 
involved are experimental and keep all the theoretical inputs fixed. The thick lines show the x^ 

^Note that this is a very conservative estimate. The total systematic error in the SK solar neutrino data sample 
is only 3.5% [19]. 
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for the case where the N"^'^'"'^ correspond to fi;(true) = 1 and hence to a "true" FD spectrum. The 
thin hnes give the corresponding for the case when K(truc) — —1 and the BE spectrum is true. 
For all the cases we generated the "data" for Tp^(true) = 5 MeV, while for Ta.(true) we assume 
four different values displayed in the figure. If we define^ our 3(j limit as given by Ax^ = 9, then 
the figure shows that for a true FD spectrum, SK could in principle restrict k > 0.2 at 3a, for 
Ti>^ (true) = 5 MeV and Ta;(true) = 8 MeV^. In particular, the false BE distribution {n = — 1) 
can be ruled out comprehensively at about 9a level. If on the other hand BE distribution was 
true, then all k, > —0.5 could be ruled out at 3a and the false FD distribution {k, — 1) could be 
disfavored at about 9a level. 

In the discussion above we had assumed that all theoretical parameters involving the prediction 
of the SN neutrino fluxes were absolutely known. However, this is far from true. In fact, the error 
coming from our lack of correct modeling for the supernova neutrino energies and luminosities is 
much larger than the experimental errors. In the last few years it has been realized [9] that the 
nucleon brcmsstrahlung NN NNv^^x as well as v^Ve '^x^x along with its inverse reaction 
are important and were not included in the earlier simulations of the neutrino transport inside 
the supernova. While the flrst process reduces the difference between and T^^, the second ones 
increase the v^^ luminosity. We note that these uncertainties could also be correlated. In what 
follows, we take into account the uncertainties on the predicted neutrino fluxes by keeping both the 
total luminosity and the temperatures Tp^ and completely unconstrained, thus allowing for 
all possible values for these parameters. This method therefore can be used to gauge the potential 
of SK to determine the true energy distribution of the SN neutrinos, irrespective of our (lack of) 
knowledge of the initial SN fluxes. Henceforth, we will distinguish between the "true" values of the 
parameters (denoted as Tp^(true) etc.) at which we generate our projected SN neutrino data set 
and the fitted values of those parameters (denoted as Tp^ etc.), obtained through the minimization 
of the function defined in Eq. (7). 

Fig. 5 shows the as a function of Tp^(true) for three assumed values of Ta.(true). The 
thick hnes in the figure show the case where the "data" N°^^'"^ is generated assuming a true FD 
distribution for ^'^(true) = 3 x 10^^ ergs and values of T'p^(true) between 4.5-7 MeV and for 
three different values of Tx{true) = 6, 7 and 8 McV shown by red dotted and green dashed and 
blue dot-dashed lines respectively. This data set is then fitted with jV*^'^°''^ calculated using the 
wrong BE distribution and with Tp^, and Eb allowed to take any possible value which gives the 
minimum x^ in the fit. The thin lines in the figure show the corresponding minimum x^ when the 
true distribution appearing in N°^^'"'^ is BE and the wrong distribution chosen in Nf^^'^^ is FD. We 
note from the figure that if FD is the true distribution, then the SN data set in SK can be used to 
rule out the possibility of neutrinos obeying the BE statistics at the 5(7 level for Tpg(true) > 4.8, 
5.0 and 5.2 MeV for Tj.(true) = 6, 7 and 8 MeV respectively. If the BE distribution for the 
neutrinos was true, then the wrong FD statistics could be ruled out at ba for Tp^(true) > 6.0, 5.9 
and 6.1 MeV, when Tj;(true) = 6, 7 and 8 MeV respectively. The wrong distribution can be ruled 
out at more than 3a for both cases for all plausible values of Tp^ . 

^Throughout this paper we will define our na limit as given by Ax^ = . Note that this is just a definition 
we assume here for the sake of comparison of the C.L. we obtain for the difi^erent types of fits. The C.L. in general 
should depend on the number of constrained degrees of freedom in the theory. 

^There is very little dependence on the value of T^. 
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The C.L. with which any given value of k can be ruled out for either the case where FD (thick 
hnes) or BE (thin lines) is the true distribution, can be seen from Fig. 6, where we have plotted 
the as a function of k. For all the curves we take Tp^(true) = 5 MeV, while for Tj.(true) we have 
assumed four values: 5 (black solid lines), 6 (red dotted lines), 7 (green dashed lines) and 8 (blue 
dot-dashed lines) MeV. The method of analysis is similar to that used for Fig. 5. In particular, 
for the case of true FD (BE) distribution shown by the thick (thin) lines, the data set N"^^'"'^ is 
generated for «;(true) = 1 (fi;(true) = —1), ^^^(true) = 3 x 10^^ and the fixed values of Tp^(true) 
and T'j.(true) mentioned above. This data set is then fitted with jv*'*^'"'^ calculated for each k, in 
the range [—1, 1], with the SN parameters allowed to vary freely. For the true FD distribution, we 
can limit allowed k > —0.39, —0.45 and —0.55 at 3a. for Tj,(true) = 6, 7 and 8 MeV respectively. 
For BE as the true distribution, k. < 0.13, 0.21 and 0.38 is allowed at 3cr, for Ta,(true) = 6, 7 and 
8 MeV respectively. 

From the Fig. 6 we note that when a theoretical distribution with a wrong k is used to fit the 
data, the obtained for the true FD distribution is in general larger than that for the true BE 
distribution. This happens because by changing the (anti) neutrino temperatures and luminosity, a 
spectrum with mixed statistics can reproduce the data set more easily when the latter corresponds 
to a BE distribution. We have checked that the spectral shape is the most important factor in 
this case. Indeed if we analyzed the data set on total (Pg +p) observed rate in SK rather than the 
positron energy spectrum, the sensitivity to k would get completely lost since the total rate can 
always be reproduced by a theory with any k when Eb, Tp^ and are allowed to vary freely. 

Fig. 7 shows the we can expect when a data set corresponding to a particular case of mixed 
neutrino spin-statistics is fitted with either a pure FD or a pure BE distribution. We generate the 
data set corresponding to each of values of — 1 < «;(true) < 1 and fit it with either a theory with 
FD distribution (thick lines) or a theory with BE distribution (thin lines). We can draw similar 
inferences from this plot as obtained from Fig. 6. 

Throughout our discussion so far we have played down the impact of the oscillation parameters. 
We generated our data set at the benchmark values of Am|i = 8 x 10~^ eV^, sin^ 9i2 — 0.3, 
A.mli = 2.1 X 10~^ eV^ and sin^^i3 = and kept them fixed in the fit. For the current allowed 
range of neutrino oscillation parameters, the survival probability for the SN (anti)neutrinos depend 
on Amli and sin^ 9i2 only for (anti) neutrinos crossing the earth. Since we assume that the position 
and time of the SN in the galactic center with respect to our detector SK will be such that neutrinos 
do not cross the earth matter, we are justified in keeping Am|^ and sin^ 9i2 fixed at any value of 
their current allowed range. However, the oscillations inside the SN matter depend crucially on 
the values of Amg^ and 613 [15]. Most importantly, for a given value of 6*13 (> 10"®), whether large 
flavor oscillations appear in the neutrino or the antincutrino channel is determined by the sign of 
Am|i, i.e. the neutrino mass hierarchy. The sign of Am^i is typically expected to be determined 
using matter effects in the 1-3 channel. For very large values of 6'i3, synergies between the T2K 
and NOz/A experiments could somewhat indicate the sgnlAmli) [20]. However, an unambiguous 
measurement would require either a beta beam facility or a neutrino factory [21]. Resonant matter 
effects in the 1-3 channel encountered by atmospheric neutrinos can be exploited to probe the 
neutrino hierarchy both in water Cerenkov and large magnetized iron calorimeter detectors [22]. 
Recently some novel ways of probing the mass hierarchy requiring very precise measurements and 
using the "interference terms" between the different oscillation frequencies have been proposed [23] . 
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Neutrinoless double beta decay experiments have the potential to provide us with the neutrino 
mass hierarchy even for vanishing ^13 [24, 25]. The uncertainty on the value/limit on ^13 is also 
expected to reduce in the future. The upper limit could be improved further to sin^ 613 ^ 0.0025 
[26] by the combination of the next generation beam experiments T2K and NOz/A, as well as by 
the second generation reactor experiments [27]. It should be borne in mind that any one of these 
experiments could even measure a nonzero 6*13, if the true value of 6*13 happens to fall within their 
range of sensitivity. 

In our analyzes so far, we had kept ^13 fixed at zero. For this case the oscillation probability is 
independent of the neutrino mass hierarchy for neutrinos not crossing the earth. Hence, our study 
so far is valid irrespective of the neutrino mass hierarchy. However, there is no reason to assume 
that the mixing angle ^13 = 0. In fact, large number of theoretically attractive models predict that 
^13 is large and possibly close to its current upper limit. If the true value of ^13 is indeed large then 
we expect big matter induced fiavor oscillations in the SN for the antineutrino (neutrino) channel, 
when the neutrino mass hierarchy is inverted (normal). For 10~^ ^ sin^ ^13 ^ 10~^ these matter 
induced conversions in the antineutrino (neutrino) channel for inverted (normal) hierarchy are also 
energy dependent, causing larger suppression for antineutrinos (neutrinos) with smaller energies 
[15]. To take into account the impact of our lack of knowledge of the true value of ^13 we present the 

as a function of k, in Fig. 8, which is similar to Fig. 6 but with ^13 allowed to vary freely as well. 
The survival probability for Pg is same for both normal and inverted hierarchy for sin^ ^13 ^ 10~^. 
For the normal hierarchy the survival probability for Pg remains the same for all ^13, while for the 
inverted hierarchy the probability reduces as sin^ ^13 increases. Therefore, by taking the inverted 
hierarchy in the fit and allowing the value of sin^ 6*13 to vary in the full range [0-0.04], we take into 
account both the uncertainty due to ^13 as well as the hierarchy [28] . A comparison of Fig. 8 with 
Fig. 6 brings out the fact that the for the case when the data set corresponds to FD spectrum 
reduces substantially when 6*13 is allowed to vary in the fit. As discussed in detail before, the SN 
neutrino spectrum depends on whether the neutrinos follow the FD, BE or mixed statistics. In 
particular, we saw in Figs. 1 and 3 that the FD distribution allows for a spectrum enriched in 
higher energy and depleted in lower energy (anti)neutrinos, compared to a spectrum with BE or 
mixed statistics. Since as noted above, for 10^® < sin^ ^^13 < 10"^ there is larger suppression in the 
lower energy end of the event spectrum due to matter induced oscillations, a theoretical BE (or 
mixed statistics) spectrum with sin^ ^13 in the range [10~^-10~'^] and with inverted hierarchy can 
be reconciled better with a data set generated for the FD distribution with ^13 = 0. This helps in 
reducing the when the FD distribution is true. However, we can see from the Fig. 8 that the 
corresponding to the true BE distribution does not change much due to the uncertainty in ^^13 and 
hierarchy. This is because oscillations further accentuate the problem if a true BE distribution is 
to be fitted with a k > — 1, except when Ti>^(true) ~ Tr(true). 

To summarize, we use the number of positron events induced in SK through the Pg capture on 
protons for a future type-II galactic supernova, to probe whether the spin-statistics is violated for 
neutrinos. We assumed that the neutrinos behave as prescribed by the standard model (including 
massive neutrinos) in all respects, except for the fact that their thermal distribution might obey 
a mixed statistical distribution rather than FD. We performed a detailed x^ analysis of the "ob- 
served" event spectrum in SK to find the allowed range of the Fermi-Bose parameter k. We have 
presented results incorporating the uncertainties stemming from the astrophysical parameters for 
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the supernova as well as those coming from neutrino oscillation parameters. 
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Figure 1: The thermal energy spectrum for the FD, MB and BE distribution of SN neutrinos with 
typical temperature of 3.5 MeV (left panel), 5 MeV (middle panel) and 8 MeV (right panel). 
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Figure 2: The neutrino spectrum expected on earth from a galactic SN for the pure FD (black 
solid line), pinched FD (red dotted line) and pure BE distribution (green dot-dashed line). The 
average energy for all the three spectrum is taken as 15.75 MeV. 
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Figure 3: The {v^+p) event spectrum expected in SK from a galactic SN. The sohd black and red 
lines show the events spectrum for FD and BE distribution respectively where neutrino oscillations 
have been taken into account. The dashed black and red lines show the corresponding spectra 
without oscillations. 
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Figure 4: The ^ function of k. The thick (thin) hues show the obtained between a 

data set for a true FD (BE) spectrum and a theoretical spectrum with mixed statistics. The four 
hues correspond to four different values of Tj.(true). For all cases we have taken a fixed value of 
Tp^{tTue) = 5 MeV. The oscillation parameters and Eb are also kept fixed. 
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Figure 5: The as a function of Tp^(true) between a data set obtained for a true FD/BE spectrum 
and a theoretical spectrum with BE/FD statistics. The thick lines are for data set with true FD 
statistics and the thin lines for data set with true BE statistics. We take three different values for 
T3;(true) and let all SN parameters to vary freely in the fit. 
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Figure 6: The ^is a function of k between a data set obtained for a true FD/BE spectrum and a 
theoretical spectrum with mixed statistics. The thick hues are for data set with true FD statistics 
and the thin hues for data set with true BE statistics. We take four different values for Tj;(true) 
and let all SN parameters to vary freely in the fit. 




K(tme) 

Figure 7: The as a function of fi;(true). Everything else is same as in Fig. 6. 
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